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We study the magnetization process of the spin-1/2 Heisenberg antiferromagnet on a 
distorted square-kagome lattice by the numerical-diagonalization method. The magnetization 
jump at one-third of the height of the saturation is examined in detail; we find that the jump 
becomes larger when a small distortion is switched on and that it is accompanied by an abrupt 
change in lines along microscopic spin directions. Our finite-size results successfully confirm 
that the magnetization jump in a spin-isotropic system is a macroscopic jump that survives in 
the thermodynamic limit and that the changes in spin directions are common to a spin-flop 
phenomenon observed in spin-anisotropic systems. 


1. Introduction 

Frustration has attracted much attention in condensed 
matter physics because it plays an important role as 
an origin of nontrivial behaviors in various systems. 
Particularly, such nontrivial behaviors are often observed 
in magnetic materials. One of them is the magnetization 
plateau. A magnetization plateau is the appearance of 
a region of a magnetic field in a magnetization process 
where the magnetization does not increase even with 
an increase in the magnetic field. This phenomenon is 
in contrast to a normal case, in which a magnetization 
process shows a smooth and signihcant increase in 
magnetization with an increase in magnetic field. Such 
a plateau originates from the existence of an energy 
gap between levels in magnetic fields owing to the 
formation of an energetically stable quantum spin state, 
for example, the spin-1 alternating spin chain^ and 
two-dimensional orthogonal dimer system.^ 

Another nontrivial behavior is the magnetization 
jump. During a jump in the magnetization process, 
the increase in magnetization is discontinuous, in 
contrast to the fact that a magnetization plateau 
corresponds to a discontinuity with respect to the 
magnetic field. As an origin of such magnetization 
jumps, a spin-flop phenomenon is well known.^ The 
phenomenon is the occurrence of an abrupt change in 
lines along microscopic spin directions while the states 
change owing to the increase in magnetic field. It is 
widely known that the phenomenon occurs when the 
system includes anisotropy in spin space.Under these 
circumstances, there was a report of an interesting 
system that shows a magnetization jump even when 
the system has no anisotropy in spin space; the system 
is the spin-1/2 Heisenberg antiferromagnet on the 
square-kagome lattice.® 

This lattice is originally introduced in the study of 
the relationship between the spin model on this lattice 
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and the eight-vertex model; ^ numerical-diagonalization 
studies of the quantum system on this lattice were 
carried out®’® although the above jump was not 
recognized within these studies. The lattice shares several 
characteristics with the kagome lattice; for example, 
the coordination number z = 4 from a vertex of each 
lattice is the same. The corner sharing of neighboring 
local triangles is also the same. We hereafter call the 
structure of each unit cell a shuriken from its shape. 
The most important difference between the kagome 
and square-kagome lattices is whether or not all the 
vertices are equivalent with respect to the polygons 
that surround a focused vertex. In the square-kagome 
lattice, vertices are divided into two groups: one that is 
a vertex of the small square inside the local structure 
of a shuriken and one that is not (see Fig. 1). Let us 
call a vertex site in the former (latter) group the a (/3) 
site hereafter. This difference induces nontrivial different 
behaviors between the magnetization processes of these 
two systems. In the kagome-lattice antiferromagnet, 
there is a region similar to a magnetization plateau 
at one-third height of the saturation. Just outside 
of this height, the change in the magnetization is 
continuous but the magnetization process shows a 
nontrivial critical behavior, which is certainly different 
from a magnetization plateau of the well-known type in 
a two-dimensional system. Although there appears 
a magnetization plateau at the same height in the 
square-kagome-lattice antiferromagnet, the plateau was 
found to be accompanied by a discontinuous change in 
the magnetization at the edge of the higher-held side in 
the magnetization processes of hnite-size clusters. The 
discontinuity strongly suggests a magnetization jump. 
The local magnetization also shows its discontinuity, 
which suggests that a behavior similar to a spin-hop 
phenomenon happens. A common behavior in hnite-size 
clusters is observed in other cases: the Cairo-pentagon 
lattice,^®’the distorted kagome lattice,^®’and the 
shuriken-honded honeycomb lattice.^® 

However, to conhrm that this discontinuous behavior 
of hnite-size clusters in the square-kagome-lattice 
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antiferromagnet corresponds to a true magnetization 
jump owing to a spin-flop phenomenon, the following 
two viewpoints should be clarified. One is whether or 
not the jump certainly survives in the thermodynamic 
limit. In a normal magnetization process without any 
anomalies, the increase in magnetization of a finite-size 
cluster must be 5M = 1; in the square-kagome-lattice 
antiferromagnet, a discontinuity of 5M = 2 was found. 
Note here that a discontinuity of 5M = 2 can appear 
in the spin-nematic state.Since the spin-nematic 
state is the two-spin magnon state, the discontinuity of 
the spin-nematic state is not a macroscopic behavior. To 
exclude the possibility that a discontinuity of 5M = 2 
in the square-kagome-lattice antiferromagnet indicates 
the spin-nematic state, it is sufficient to observe a 
discontinuity of 5M > 3. The other point is how 
the transverse spin component behaves. The local 
magnetization examined in Ref. 6 provides information 
on a component along the external field, namely, the 
z-axis. If the transverse spin component between both 
sides at the observed jump shows a significant difference, 
the relationship between the magnetization jump and 
the spin-flop phenomenon is established. Under these 
circumstances, the purpose of this study is to clarify 
these two unresolved issues in the square-kagome-lattice 
antiferromagnet by taking distortion into account 
additionally to the undistorted square-kagome-lattice 
case. 

Candidate square-kagome-lattice antiferromagnet 
materials have not been reported so far to the best 
of our knowledge. On the other hand, the discoveries 
of candidate kagome-lattice antiferromagnet materials 
such as herbertsmithite,^®’volborthite,^°“^^ and 
vesignieite,^®’^"^ have accelerated experimental studies 
as well as theoretical ones,^®“^^ leading to a deeper 
understanding of the physics of the kagome-lattice 
antiferromagnet. From these studies, it is gradually 
clarified that the materials have significant deviations 
from the ideal situation of the kagome-lattice Heisenberg 
antiferromagnet. In volborthite, for example, magnetic 
interactions deviate from the uniform situation in the 
kagome lattice because of the orbital directions despite 
the fact that atomic positions are located at the vertices 
of the kagome lattice. To understand the candidate 
materials well, therefore, effects due to the deviation 
from the ideal situation are worth studying, which is a 
background of this study. 

This paper is organized as follows. In the next section, 
the model that we study here is introduced. The method 
and analysis procedure are also explained. The third 
section is devoted to the presentation and discussion 
of our results. We first discuss the change in local 
magnetizations of the state at one-third of the height 
of the saturation. The discussion will clarify which phase 
the undistorted square-kagome-lattice case is located in. 
After that, the magnetization process of larger clusters is 
presented. The spin-spin correlations are also examined. 
In the final section, we present our conclusion. 



(d) (e) 




Fig. 1. (Color) Structure of the square-kagome lattice. In panel 
(a), a. and ^ sites are illustrated; the broken-line square 
represents a unit cell of this lattice. The shapes of finite-size 
clusters studied in this work are shown for (b) Ns = 24, (c) 
Ns = 30, (d) Ns = 36, and (e) Ns = 42. 

2. Model Hamiltonians and Method 

The Hamiltonian that we study in this work is given 
by H = Ho + Hzeeman, where 

Ho= ^ JiS,-Sj+ Y. (1) 

{i,j) ,i^OL,jGl3 {i,j),i^OL,jGcx 

and 

Hzeeinan = —h 'Y, ■ (2) 

3 

Here, Si denotes the S' = 1/2 spin operator at site *, 
where the sites are the vertices of the square-kagome 
lattice shown in Fig. 1. The spin operator satisfies 
= S(S + 1). The sum of Ho runs over all the 
nearest-neighbor pairs in the square-kagome lattice. 
Energies are measured in units of Ji; hereafter, we set 
Ji = 1. The number of spin sites is denoted by Ag, where 
As/6 is an integer. We treat As = 24, 30, 36, and 42 in 
this study. We impose the periodic boundary conditions 
for clusters with site As; the shapes of the clusters are 
shown in Fig. 1. Note here that the shapes for As = 24 
and 30 are squares although the latter is slightly tilted. 
On the other hand, the shapes for As = 36 and 42 are not 
squares. However, examination of systems with a large 
As would give a higher resolution of the magnetization 
process and contribute much to our understanding of this 
system. Note here that the same model including the 
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control of J 2 was previously examined in Refs. 46 and 
47. Reference 46 treated this model within an effective 
Hamiltonian based on the localized-magnon picture. 
Reference 47 treated only Ns = 24 and 30. Both of 
these studies did not examine the magnetization jump 
at one-third of the height of the saturation on which we 
focus our attention in this study. 

We calculate the lowest energy of T-Lq in the subspace 
belonging to = M hy numerical diagonalizations 

based on the Lanczos algorithm and/or Householder 
algorithm. The energy is denoted by E{Ns,M), where 
M takes an integer value down from Miow (=0) up 
to the saturation value Ms (= SNs). We use the 
normalized magnetization m = M/Ms- Some of the 
Lanczos diagonalizations were carried out using the 
MPI-parallelized code, which was originally developed 
in the study of Haldane gaps."^® The usefulness of 
our program was previously confirmed in large-scale 
parallelized calculations.Note here that the 
calculations in cases of a small M for Ns = 42 require 
the use of the K computer (Kobe, Japan) to obtain the 
magnetization process of the entire range. 

The magnetization process for a finite-size system is 
obtained by considering the magnetization increase from 
M to M -I- 1 in the held, 

H = E{Ns,M + l)-E{Ns,M), (3) 

under the condition that the lowest-energy state with the 
magnetization M and that with the magnetization M + 1 
become the ground state in specihc magnetic helds. It 
is often the case that the lowest-energy state with the 
magnetization M does not become the ground state in 
any held. In this case, the magnetization process around 
the magnetization M is determined by the Maxwell 
construction.^’^ 

We evaluate the average local magnetization dehned 
as 

= (4) 

where ^ takes a and /3. N(^ denotes the number of ^ sites; 
namely, Na = 2iVs/3 and Np = Ns/3. Here, the symbol 
(O) represents the expectation value of the operator 
O with respect to the lowest-energy state within the 
subspace characterized by a hxed magnetization M that 
we focus our attention on. Averaging over ^ is carried 
out when the ground-state level is degenerate. Note that, 
when the ground state belonging to M is not degenerate, 
the results do not change regardless of the presence or 
absence of this average. We also measure the correlation 
functions (SjS/.) and to capture well the features 

of the wave functions that are numerically obtained. 

3. Results and Discussion 

First, we present our results of the J 2 dependence of 
the local magnetization of the m = 1/3 state, which 
are shown in Fig. 2. A similar examination concerning 
the parameter dependence was carried out in cases of 
the Cairo-pentagon lattice^® and the distorted kagome 
lattice,^® and the Lieb lattice accompanied by frustrating 
interaction.One finds in Fig. 2(a) a large decrease in 



Fig. 2. (Color) Dependence of the local magnetization at m = 
1/3 on the ratio of J 2 I Panel (b) is a zoom-in view of panel 
(a) to observe the range near the critical point in detail. Black 
closed triangles and inversed triangles denote = 24 data of 
a and /3, respectively. Green open squares and diamonds denote 
Na = 30 data of a and /5, respectively. 

m£]y[ and a large increase in near J 2 /J 1 ^ 0.96. 
This marked change suggests that a phase transition 
occurs between the ferrimagnetic state for J 2 /J 1 ^ 0.96 
and the state composed of singlet spins and almost fully 
polarized spins for J 2 /J 1 0.96. The phase transition 

has already been reported in Ref. 47, which concluded 
that this is a first-order transition. To detect the behavior 
of mLM near J 2 /J 1 ^^0.96, we observe a detailed change 
in the zoom-in view in Fig. 2(b). Our results show a 
continuous change in for both Ns = 24 and 30 

although the change for Ng = 30 is much more rapid than 
that for Ns = 24. From the observation of this continuous 
behavior, the conclusion of the first-order transition may 
be premature. Note also that the continuous change for 
Ns = 30 is in contrast to the discontinuous change for 
the same Ng in the case of the Cairo-pentagon lattice^® 
despite the fact that the present model and the model on 
the Cairo-pentagon lattice share the same number of spin 
sites in each unit cell. It is difficult within studies based 
on finite-size examinations, but it should be resolved in 
future studies to determine whether the transition is of 
the second or first order. An important point clarified 
in this study is that J 2 /J 1 = I is different from the 
transition point 0.96 and that J 2 /Ji = 1 is in the state 
composed of singlet spins in the a sites and almost fully 
polarized spins in the /3 sites. Hereafter, in the present 
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Fig. 3. (Color) Magnetization process for J 2 IJ 1 = 1-02. Black 
triangles, green diamonds, blue squares, and red circles denotes 
data for Ns = 24, 30, 36, and 42, respectively. 



Fig. 4. Average local magnetization at m = 1.02. Pluses and 
crosses denote Ns = 24 data for a and 13^ respectively. Diamonds 
and squares denote Ns = 30 data for a and respectively. All 
the states are realized for Ns = 24, but not for Ns = 30. For 
Ns = 30, closed symbols represent data for the stably realized 
states, while open symbols denote data for the unstable states 
during the magnetization jump. 

study, we focus on the case of the side of J2/J1 that is 
greater than the transition point J 2 I J\ ~ 0.96. 

Next, we examine the magnetization process in the 
phase of the J2/J1 = 1 side. Our results are shown 
in Fig. 3 for J2/J1 = 1.02. This parameter is chosen 
within the same side as J2/T1 = 1 as a value close to 
J2I Ji = 1- Some features are shared between J2I J\ = 1 
and 1.02; a detailed comparison between the results 
of J2/J1 = 1 and 1.02 will be carried out later. The 
magnetization processes in Fig. 3 show a clear existence 
of the TO = 1/3 plateau. No jump appears for Ns = 24 at 
the edge of the higher-held side of the m = \l'i height. 
For Ns > 30, on the other hand, magnetization jumps 
are clearly observed. The appearance of such jumps is 
independent of iVs, at least within Ns > 30 of the 
present hnite-size-system study. In Fig. 4, the average 
local magnetizations for the square-cluster systems are 
shown. One easily hnds that discontinuous behaviors 



h/J 



h/J 

Fig. 5. (Color) Zoom-in views of the magnetization processes for 
(a) J 2 /J 1 = 1.02 and (b) J 2 /J 1 = 1. Symbols are the same as 
those in Fig. 3. The broken lines represent the results before the 
Maxwell construction is carried out. 

in a and /? sites appear between M = (1/3)Mg and 
M > (l/3)Ms. At M = (l/3)Ms, our results suggest that 
a spin at a /3 site reveals almost a full moment along the 
z-axis and that spins at a sites form a singlet state in 
the small square inside each shuriken. 

Let us next compare between the results of the 
magnetization processes for J 2 /J 1 = 1 and 1.02. 

Figures 5(a) and 5(b) show zoom-in views of the 
magnetization processes around the edge of the 
higher-field side of the to = 1/3 height. The most marked 
difference between the two cases is the skip SM in the 
results for Ns = 42; we observe 5M = 2 for J 2 /J 1 = 1 
whereas we observe 6M = 3 for J 2 /J 1 =1.02. The result 
of 6M = 3 for J 2 /J 1 =1.02 is important from the 
viewpoint of the origin of the jump. A finite-size jump 
with 5M = 2 can appear in the nematic phase, 
where the two-magnon bound state is realized. The 
nematic phase is gapless so that the magnetization 
process in the thermodynamic limit is continuous. Thus, 
the finite-size jump owing to the nematic phase is only 
a finite-size effect and does not become a macroscopic 
jump in the thermodynamic limit. The present result of 
SM = 3 for Ns = 42 and J 2 /J 1 =1.02 clearly excludes 
the possibility for the appearance of the finite-size jump 
owing to the nematic phase. The skip SM = 3 also 
suggests that the jump becomes a macroscopic one in 
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the thermodynamic limit. 

Let us then estimate the skip in the thermodynamic 
limit from our finite-size results. The resolution of the 
normalized magnetization process is given by 2/Ns. If 
a finite-size system reveals a jump with SM = 2, the 
skip Sm of the normalized magnetization jump in the 
thermodynamic limit should satisfy 

2/Ns <Sm< 2{2/Ns), (5) 

while the skip is given by 

Sm > 212 /Ns), (6) 

if a finite-size system reveals a jump with SM = 3. For 
J2/J1 =1.02, we have SM = 2 for Ns = 36 and SM = 3 
for Ns = 42; Eqs. (5) and (6) suggest 

0.095 ^ (5m ^ 0.111. (7) 

For J2/J1 =1, on the other hand, we have SM = 2 for 
Ns = 36 and 42 and no jumps for Ns < 30, which suggest 

Sm 0.095. (8) 


behavior that we have observed in the m^^. Thus, it 
is confirmed that, at m = 1/3, the system forms the 
state in which the spins at /3 sites have almost a full 
moment along the z-axis and where spins at a sites form 
a singlet state in the small square inside each shuriken. 
At M = (l/3)Ms -I- 2, on the other hand, the spin state 
is markedly different. {S/Sp for i,j £ (3 is markedly 
smaller than that at M = (l/3)Ms while {S/Sp for 
i,j £a is larger. In addition, both {SfSp for i,j £ a and 
{SfSp for i, j £ j3 are significantly nonzero. The results 
concerning the transverse component strongly suggest 
that a and P spins take directions that have nonzero 
angles from the z-axis. This picture concerning the state 
of M = (l/3)Ms-|-2 is consistent with the speculation for 
J2/ Ji = 1 presented in Ref. 6 as the umbrella-type spin 
state. Finally, note that the behavior at M = (l/3)Ms-|-l 
is qualitatively the same as that at M = (l/3)Ms -I- 2, 
although it is unclear whether the m = 0.4167 state is 
stable or unstable. 

4. Conclusion and Remarks 


Therefore, one finds that the skip of the jump becomes 
larger when J 2 /Ji deviates from the critical point ^ 0.96 
clarified in Fig. 2. 


Table I. Correlation functions for Ns = 24 cluster at J 2 IJ 1 = 
1.02. We select a pair of sites i and j so that the distance between 
these sites is largest in this cluster. 




M = 

M= iMg + 1 

M = ^Ms + 2 

i,j £ a 

(SfSp 

(SfSp 

0.0005068 

0.0000012 

0.0571670 

0.0048659 

0.0498709 

0.0043312 

£ 0 

(Sfsp 

(SfSp 

0.2071796 

0.0000184 

0.0335379 

0.0192386 

0.0960927 

0.0167820 


Next, let us consider the spin direction from the data 
of correlation functions. We measure {SfSp and {SfSp 
between the pair of sites i and j, where i and j in the 
same group are the most distant in a finite-size cluster 
forming a square. The same analysis was carried out in 
Ref. 15. Within this study, the Ns = 24 cluster is the case 
that should be measured while the others are not because 
the Ns = 42 and 36 clusters do not form a square and 
because i and j as the most distant sites for the Ns = 30 
cluster are not in the same group. Our present results 
are summarized for J 2 /J 1 = 1.02 in Table 1. Note here 
that, for Ns = 24, M = (l/3)Ms, M = (l/3)Ms -|-1, and 
M = (l/3)Ms -I- 2 correspond to m = 1/3, m = 0.4167, 
and m = 0.5, respectively. From Eq. (7), it is clear that 
m = 0.5 is outside the jump, although it is unclear 
whether m = 0.4167 is outside or inside the jump because 
m = 0.4167 is very close to the upper edge of the jump. 
Therefore, it is reasonable to consider that the M = 
(l/3)Ms -I- 2 state is appropriate for observing a state 
under a field higher than the jump; on the other hand, the 
M = (1/3) Mg state is appropriate for observing a state 
under a field lower than the jump. At M = (1/3)Mg, 
only {SfSp for i,j £ P is dominant; the other quantities 
are very small. These results are in agreement with the 


In this paper, we study the spin-1/2 Heisenberg 
antiferromagnet on a square-kagome lattice with 
distortion by numerical-diagonalization calculations. We 
have found that the magnetization jump occurs at 
the higher-held edge of one-third of the height of 
the saturation as a frustration effect. An essential 
improvement from the results in Ref. 6 is that our 
large-scale calculations detect the behavior of the growth 
of the jump when the interactions in the small square 
in a shuriken become larger. Our hnite-size results 
have clarihed that a possibility of a hnite-size jump 
owing to the spin-nematic state is excluded. Our 
analysis of spin-spin correlation functions is another 
essential improvement suggesting that, in the higher-held 
region, each spin reveals its own signihcant transverse 
component, while no transverse components are observed 
in the state at one-third of the height of the saturation. 
The changes in spin directions in the spin-isotropic 
system share the behavior in the spin-hop phenomenon 
observed in spin-anisotropic systems. Magnetization 
jumps in spin-isotropic systems were reported in several 
studies.The comparison between these cases and 
the present case will also be useful for understanding well 
the occurrence of these jumps. The present consequence 
of the magnetization jump will become a clue for 
future investigations, which would contribute much to 
our understanding of frustration effects in magnetic 
materials. 
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